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1 Introduction 



Dynamical triangulation is a discrete approach to the path integral for euclidean 
quantum gravity in which the euclidean spacetimes are constructed by 'glueing' 
together geometric objects such as equilateral hypercubes or simplices. Hyper- 
cubes were used at first ^ whereas the modern formulation which developed 
independently out 2D gravity uses simplices. Results of numerical simulations of 
the commonly used 4D model were first presented in ^ . As a model of gravity 
dynamical triangulation ought to have a scaling regime where it corresponds to 
semiclassical Einstein gravity. One wishes to recover classical euclidean space- 
times, Newton's potential and the formation of gravitationally bound states. 

It is not clear that a purely euclidean formulation should be able to con- 
tain semiclassical gravity, because of the well known divergence related to the 
unboundedness of the euclidean version of the Einstein-Hilbert action. In a semi- 
classical evaluation of the euclidean path integral this can be dealt with by de- 
forming the integration over the conformal mode into the imaginary direction of 
the complex plane [Q. 

The dynamical triangulation formulation is completely regular from the start. 
The model has a phase transition as a function of a parameter which is propor- 
tional to the inverse bare Newton constant: K2 oc Gq^. The transition separates a 
phase with crumpled spacetimes and very high effective dimensionalities, from an 
'elongated' phase with effectively two-dimensional spacetimes with characteristics 
of a branched polymer p, |^ . Near the transition the model appears to produce 
classical S'^-like spacetimes, in an intermediate distance regime, and there is evi- 
dence for scaling, suggesting continuum behaviour[][^ ^. In the elongated phase 
the scaling degenerates into a branched polymer version [^]. 

It is of great interest to find out if this scaling region can be described by an 
effective action of the type 

S..^/</'xV5(^-^ + Cfi^ + -). (1) 

The branched polymer phase is presumably an expression of the conformal mode 
instability of the —R/G-pait in this effective action at scales of order while 
the crumpled phase may correspond to negative G. 

Recent evidence suggests that the transition is of first order [0, ||, instead 
of second order as thought previously. The first order nature of the transition 
need not stand in the way of continuum behaviour as the examples of gauge- 
Higgs models show. These models have a first order phase transition also in 
the continuum treatment. For dynamical triangulation we suggested previously 

^We use the phrase 'continuum behaviour' rather than 'continuum hmit' to allow for situa- 
tions in which the lattice distance is very much shorter than any physical distance, but for which 
the continuum limit would imply unwanted features, such as noninteraction due to triviality (a 
typical example is the 4D Ising model formulation of theory). 
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that continuum behavior may be automatic, as in two dimensions^]. It is a 
priori possible that continuum analysis of the effective theory (0) also leads to a 
first order transition between phases with positive^ and negative 1/G. The above 
interpretations are very speculative and more analytical as well as numerical work 
has to be done to establish or reject it. 

It would be very interesting if we could measure the attraction of two sources 
at a fixed distance. This has been pursued in the Regge calculus formulation 
of simplicial quantum gravity [10|. Such measurements could then be compared 
with the simple Newtonian law or with quantum corrections to this law, which 
have for instance been calculated in and [O. In dynamical triangulation 



a computation of the gravitational potential appears nontrivial, because it is 
difficult in a fluctuating spacetime to keep two heavy test masses at a fixed 
distance. 

However, the formation of bound states out of two test particles can be com- 
puted in a way that is customary in lattice field theory and we shall report on 
such computations in this paper. Preliminary results have appeared earlier in 
[|l^ and We see this work as an important ingredient for the interpretation 
of 4D dynamical triangulation. 



2 Binding in the continuum 

Consider a free scalar field with bare mass itlq in a quantum gravity back- 
ground. In continuum language, the euclidean action of this system is a sum of 
a gravitational and a matter part 

S = S[g] + S[g,<P], (2) 

^[^] = T7^ [d'x^{2Ao-R), (3) 
IbTrGo J 

S[g,<P] = j d'x^(^-g^''d,<pd,<P + \mli^, (4) 

where Aq is the bare cosmological constant, R is the scalar curvature and Gq is 
the bare Newton constant. 

We take as a test particle here, i.e. the back reaction of the field (f) on the 
metric is not taken into account. In lattice QCD this approximation is often 

^We recall also an analogy with Z{n) gauge theories, which have a Coulomb phase for 
n ^ 5, with massless photons in a region of parameter space, without the need for tuning to 
a critical point. Gauge-Higgs models of this type can have furthermore a first order Coulomb- 
Higgs transition The models approach the U(l)-Higgs model as n — > oo. Because of 
triviality a continuum limit will be noninteracting, but the models may still be regarded (for 
finite sufficiently large n) as a nonperturbative formulation of the abelian Higgs model in the 
continuum. 

■^In dynamical triangulation the transition actually appears to occur at small positive l/G, 
since the spacetimes found in its neighbourhood are 5'*-like, cf. sect. 5. 
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called the quenched approximation, or valence quark approximation, because 
it neglects diagrams with internal quark (in our case 0) loops. For not too 
light quark masses it turns out to give good results (for a discussion see e.g. 
|15[]). In dynamical triangulation the inclusion of a scalar field is no problem in 
principle (and appears to have little influence on the gravity sector of the theory) 
but the enlargement of parameter space by mo is computationally costly. 
A continuum calculation of the gravitational attraction of a scalar field in the 



quenched approximation was carried out in |]TT . 

We will use the following notation for expectation values of an observable A. 
On a fixed background geometry we can average over configurations of the matter 
field 

/ V<f)Aexp{-S[g,<j)]) 
(A), = ^—r , (5) 



V(j) exp(-S'[5(, 



and we can average over metrics 



VgAeM-S[g]) 

{A), = , (6) 

Vg exp(-5[^]) 



The quenched expectation value is then 

{A) = m^),. (7) 

We next consider propagators. In a fixed geometry, the one particle propaga- 
tor, denoted by G{x, y; g), is defined as 

Gix,y;g) = {(f)^(f)y)^, (8) 

and the two-particle propagator is simply the square of the one particle prop- 
agator. Letting the metric fluctuate, we take the average of the propagators 
over the different metrics. Because of reparametrisation invariance, the average 
{G{x, y; g))g can only depend on whether x and y coincide or not. Therefore, we 
look at averages at flxed geodesic distance r. 



d x^/g G{x, y; g)5{d{x, y) - r) 

G{r) = ( ^ -r ) , (9) 

d'^Xy/g6{d{x,y) - r) 



where d{x, y) is the minimal geodesic distance between x and y. By translation 
invariance, the resulting G{r) is independent of y. Similarly, we can deflne the 
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geometry average of the two-particle propagator as 



G(^) (r) = ( ^ , ) . (10) 

d^Xy/g6{d{x,y) - r) 



In (§) we first averaged over the volume at distance r from the point y and 
then averaged over the metrics. Alternatively we can integrate G{x, y; g) over 
the volume at distance r, average over metrics, and then divide by the same with 
G ^ 1, 



d G{x, y\ g)S{d{x, y) - r) 

G{r) = ^ ^ (11) 

d'^Xy/g6{d{x,y) - r) 



A similar ambiguity arises with purely geometric correlators IjT^- Our physical 
intuition tends to favour the form (P), and we shall later use its analogue in 
dynamical triangulation. However, if there is no pair x, y with d{x, y) = r for a 
given metric g the expression becomes mathematically ill defined, whereas (pi]) 
has no such problem. For the pure geometry correlators this does not happen in 
practise for reasonable r and the difference between (^) and (^) appears to be 
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For a massive particle, we expect the propagator (^ to fall off exponentially 

as 

G(r) = Zr"exp(-mr), (12) 

with some power a and renormalised mass m, which in general will not equal the 
bare mass tjiq. This expression neglects finite size effects and should be modified 
when looking at distances comparable to a typical length scale in the system. 
The two-particle propagator will behave similarly as 

G'(2)(r) = Z(2)r^exp(-Mr), (13) 

where M is the energy of the two-particle state. If this energy turns out to be less 
than two times the mass of a single particle, the difference can be interpreted as 
a binding energy between the particles. This would show gravitational attraction 
between them. 



3 Propagating in constant curvature 

Since the spacetimes to be used in the binding energy computations have topology 
5"^ and the average spacetime is expected to be homogeneous, it will be useful to 
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know the properties of scalar field propagators on spaces of constant curvature. 
To calculate such propagators we have to solve the equation 

(□-m2)G = 0, (14) 

with the boundary condition G (47r^r^)~^ as the geodesic distance r ^ 0. We 
assume spherical symmetry, i.e. G depends only on r. We need to distinguish 
three cases: positive, zero and negative curvature. For a space with constant 
positive curvature, i.e. a four-sphere, equation (|T^) can be written as 

(9^ + 3 cot x)d^G - m^rlG = 0, (15) 

where x = r/ro, with tq the curvature radius. The substitution z = cos^(a;/2), 
which was used in [|19[, turns this into a hypergeometric equation. 



[^(1 - z)dl + (2 - 4:z)d, - mV^] G{z) = 0. (16) 
For numerical evaluation we found it easiest to use a series representation 

oo 

G(^) = iV^Cfc^^ (17) 

fc=0 

Setting Co = 1, the differential equation fixes the Ck as 

(A; + 2)(fc-l)+mVo\ 

= k{kTY) ''-'^ ^^^^ 

If we demand that the singularity at the origin goes like (47r^r^)^^, this fixes the 
normalisation to be 

N= / , (19) 

WttTq cosh ( vr Wm^rQ — | j 

One can easily check that this function also conforms to the differential equation 
at the point opposite the origin, where x = it. When the mass is so small that 
m^rQ < 9/4, the above formulas may be analytically continued such that 

2^2 o 



m r, 



N = . " . • (20) 

IGitTq cos [tc\/^ — rn'^rQ] 

The negative curvature case gives similarly 

{d^ + 3 coth x)d^G -m^lG = 0. (21) 
Using the substitution z = sechx, we get 

[(1 - z^)z^dl + (-2 - 2z^)zd, - m^l] G{z) = 0. (22) 
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In this case we use the series 

oo 

G{z)=NY,CkZ^^^, (23) 



A;=0 

and the recurrence equation is now 



^ (fc + /x-2)(A: + /i-l) ^24) 
(A; + /i)(fc + /i — 3) — rn?rl 



We can read off that we have to choose 



^^=\ + \j\ + m^rl (25) 

to avoid generating negative powers of z. Setting again cq = 1 and demanding 
that G goes like (47r^r^)~^ for small r gives us the normalisation 

r(/i/2)r(/x/2 + 1/2) 

= 4vrV§r(;.-l/2) • 
The fiat case is the easiest. Using the spherical symmetry results in 

dr + - ) drG{r) - m^G{r) = 0, (27) 



which can be directly solved in terms of a Bessel function as 

^(^) = (28) 

using the condition that G{r) goes like (47r^r^)~^ for small r. 

We have plotted an example of the three cases in figure |I|, using a curvature 
radius tq of 13 and a mass m of 0.15. These parameters where chosen for later 
comparison with dynamical triangulation results. 



4 Dynamical triangulation 

In the dynamical triangulation model of four dimensional euclidean quantum 
gravity the path integral over metrics on a certain manifold is defined by a 
weighted sum over all ways to glue four-simplices together at the faces p|, |^. 
This idea was first formulated in [Q], using hypercubes instead of simplices. In 
four dimensions the analogue of the continuum gravitational action is 

S[9] = d'x ^ {2A, - R) (29) 

^ KiNi-K2N2, (30) 
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Figure 1: Propagators on spaces with constant curvature R = 12/rQ (upper), 
R = (middle) and R = —12/rQ (lower), with tq = 13 and m = 0.15. 



where N2 and are the number of triangles and four-simplices respectively. 

The partition function of the model for a fixed volume (fixed number of four- 
simplices A^4) is given by 

Z{N,K2)= Yl exp(/€2iV2). (31) 

T{N4=N) 

The sum is over all ways to glue N four-simplices together, such that the result- 
ing complex satisfies the manifold condition, with some fixed topology which is 
usually (as well as in this article) taken to be S^. The coupling constant k,2 is 
proportional to the inverse of the bare Newton constant: 

K2 = V2/4G0, (32) 

where V2 is the area of a triangle. 

It turns out that the model has two phases |2y, ^ |2^, For low ^2 the 
system is in a crumpled phase, where the average number of simplices around a 
vertex is large and the average distance between two simplices is small. In this 
phase the volume within a distance r appears to increase exponentially with r, 
a behaviour like that of a space with constant negative curvature. At high K2 
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the system is in an elongated phase and resembles a branched polymer. As is 
the case with a branched polymer, the (large scale) internal fractal dimension 
is 2. The transition between the two phases occurs at a critical value K2 which 
depends somewhat on N . The phase transition appears to be of first order |^. 
At the transition the spaces behave on the average in several respects like the 
four dimensional sphere 0. Some of the evidence for this will be reviewed in the 
next section. 



5 Spacetimes near the transition 

We have performed numerical simulations of four dimensional dynamical trian- 
gulation, according to the partition function (0). We used systems of about 
32000 simplices and the topology of the four-sphere. To keep the number of sim- 
phces around the desired value, we added a quadratic term to the action as was 



described in 20, 23, 24 



In ref. we studied the euclidean spacetimes generated in numerical sim- 
ulations of the model by measuring the number of simplices N'{r) at geodesic 
distance r. (The geodesic distance d^y between two simplices with centres x and 
y is defined as the minimum number of links on the dual lattice between x and 
y.) We fitted N'{r) in an intermediate distance regime by a form csin'^~^(r/ro) 
corresponding to a (i-sphere of radius tq. For K2 near the transition i^2{^) ^^is 
gave (i ~ 4, which we took as evidence for classical S*^ behaviour at the distance 
scale involved. Below the transition (i.e. ^2 < 1^2) d rises steeply to large values 
while above the transition d falls rapidly to the branched polymer value d = 2. 

To get a feeling for the geodesic distances which will appear later in binding 
energy computations, we show in fig. |] the quantity N'{r) for a system of = 
32000 simplices at K2 = 1.255 and 1.259. The value 1.255 is very near K2(A^) as 
defined by the position of the maximum in the susceptibility 



dHnZ{N,K2) ^ , 



No^\ /N, 



(33) 



The sin'^"^(r/ro)-fit gives d = 4.24(4) and is seen to be reasonable in the region 
6 ^ r ^ 24. The value 1.259 is our k,2 closest to the phase transition for which the 
data fit the four-sphere also reasonably well with d = 3.67(7), in the region 5 ^ 
r ^ 25. The overall shape of N'{r) is quite asymmetrical, which is presumably 
due to branching fluctuations at larger distances. 

In we also introduced an effective curvature -Reff('^) to describe the curva- 
ture at scales much larger than the lattice scale. The argument r of Res{f) is 
meant to approach zero provided that, and as long as, -Reff('") is stationary. This 
does not seem to happen in the elongated phase, while near the transition and in 
the crumpled phase we found that -Refr had indeed a stationary point (minimum). 
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Figure 2: Fits to iV'(r) of the form a sin'' ^(r/ro) near the transition. Upper curve: 
K2 = 1.255, d = 4.24(4), vq = 13.35(7); lower curve: K2 = 1.259, d = 3.67(7), 
ro = 14.2(2). 

For smaller r-values Res rises steeply when r gets smaller. We called this region 
of large Res the 'planckian regime'. 

Fig. ^ shows Res for N = 32000 and several values of K2. For k,2 < 1-255 
we see negative curvature at the minimum of Res- A homogeneous space with 
constant negative curvature is unbounded, so the maximum in N'{r) is evidently 
due to finite size effects. We take the 'planckian regime' as the small r region, 
roughly ending at the minimum of i?efr- For k,2 = 1-255 the effective curvature 
is nearly zero at the minimum, while positive for larger K2. The minimum has 
roughly turned into an inflection point at K2 = 1.259, and it has disappeared 
altogether for larger K2- 

Fig. suggests positive curvature with a curvature radius of tq ~ 21/(7r/2) ^ 
13. This corresponds to some average near the minima in fig. |^. The minima 
will show somewhat smaller curvature: -Reff('"min) ~ 0, 0.045, for K2 = 1.255, 
1.259, respectively. The latter minimum corresponds to a curvature radius ro = 
1^12/0.045 ~ 16, which is not unreasonable compared to the previous 13. 

We conclude that for ^2 = 1.255 — 1.259 the spacetimes are on the average 
near S"^ in the distance regime 6 ^ r ^ 24. 
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Figure 3: Effective curvature plots for k.2 = 1.250-1.260. 



6 Binding in dynamical triangulation 



On each dynamical triangulation configuration we calculated the propagator of 
the scalar field 

G..„ = i-D + ml)-], (34) 



'xy \ >— ' I '"'O/xj/) 

using the algebraic multigrid routine AMG1R5, where x is an arbitrary origin. 
The discrete Laplacian is defined as 



(□) 



xy 



1 ii X and y are nearest neighbours, 
—5 if X = y, 

otherwise. 



(35) 



The 5 in the second line arises as the coordination number of a four-simplex, i.e. 
a four-simplex has five neighbours. 

We can then calculate G{r) and G^'^\r) by averaging Gxy respectively its 
square over all points y at distance r from the origin x, and then over origins and 
configurations. 



G{r) = 



(36) 
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(37) 



(we have not indicated the average over origins x). Notice that (^) corresponds 
to (P) in the continuum. 

To improve the calculation of the binding energy, we can try to use what are 
called "smeared sources" . The use of these smeared sources can improve the data 
by increasing the contribution of the ground state and decreasing the contribution 
of the excited states. Instead of using {^7\j , we can calculate G{r) by averaging 



Gxy over all points at distance r from the origin x, and only after taking this 
average, square it for the calculation of G^'^\r): 

X/y ^xySdxy,r 




(38) 



This corresponds to taking the propagator from a source that is not a single 
point, but a complete shell around the origin. Such a source may have a bigger 
overlap with the ground state wave function and a smaller one with the excited 
wave functions. For a discussion of the use of these smeared sources in QCD, see 
e.g. 
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Because the average of the square of a fluctuating quantity is always greater 
than the square of its average, it is obvious that G^'^\r) > G{r). This does not 
yet imply anything about the way they fall off. In particular it is not guaranteed 
that the energy of the two-particle state is less than twice the energy of the one 
particle states. 

In figure |^ we see the results for four different bare masses. Each pair of lines 
corresponds to one bare mass. In each pair the upper line is G^'^\r) (using (|37D ) 
and the lower line is G'(r)^ of eq. (|36|) . We used 144 configurations recorded every 
5000 sweeps (1 sweep = N accepted moves). For each of the masses we used 120 
origins per configuration. The coupling constant K2 = 1.255, which is the lower 
of the two K2 values used in the N'{r) figure |^. 

There is clearly a difference in slope between the lines in each pair. This 
shows that the energy of the two particle state is less than two times the mass of 
a single particle and consequently that there is a positive binding energy between 
the particles. 

Using this data we can measure the renormalised mass m, by assuming a long 
distance behaviour of r"exp(— mr). The results are (in parenthesis is the value 
of K2) 



mo 


m(1.255) 


m(1.259) 


0.0316 


0.14 


0.12 


0.1 


0.29 


0.27 


0.316 


0.60 


0.58 


1 


1.21 


1.20 



(39) 
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Figure 4: The two-particle propagator and the square of the one particle propa- 
gator versus the geodesic distance for four different bare masses mo- The vertical 
scale is logarithmic. k,2 = 1-255, which is very close to the transition. 



It was argued in |]20| that the physical mass should vanish at zero bare mass 
and that therefore the renormalisation would be only multiplicative. Our data 
seem to show that the relation is more complicated. Increasing mg by a factor of 
\/T0 «i 3.16 increases m by a factor of about 2.1. 

Comparing figure |l| with figure we see that the long distance behaviour is 
indeed similar, being an exponential. The finite size effect of the S'^-like curvature 
is apparently negligible, except perhaps for the smallest mass. The short distance 
behaviour is quite different. The propagators in figure ^ curve downward towards 
the origin, while the free propagators shown in figure ^ curve upwards due to the 
1/r^ behaviour. The curving downward is unusual, because a propagator is inter- 
preted as a sum of decaying exponentials corresponding to the ground state and 
various excited states. A closer look shows that the downward curvature occurs 
for distances smaller than about r = 5, which is roughly the end of the 'planckian 
regime' mentioned in the previous section. As we have seen in figure |^, larger 
positive curvature means that the propagator decreases more slowly. Therefore, 
such a planckian regime may cause the propagator to decrease less at smaller dis- 
tances where the effective curvature is large than at the longer distances where 
the effective curvature is small. 
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Figure 5: The effective binding energy Eb as a function of tlie geodesic distance 
for four different bare masses at K2 = 1.255. Tlie lines are fits tot tlie data witfi 
tlie form (H^). 



Using tliese data, we can now estimate tlie binding energy of the particles. 
From (|T^) and (|TB|) we have 



Eh = 2m-M, r ^00. (41) 

As we cannot use infinite distances, we will consider the effective binding energy 
Ei,{r) for finite r and look whether this expression becomes constant. 

Figure |^ shows this quantity as a function of the geodesic distance, using 
the smeared estimator ( p8D for G^^\r). The four curves again correspond to the 
four different bare masses in figure |. To avoid the correlations between origins 
on the same configuration and between points at the same distance of such an 
origin influencing the error bars, we first averaged all the measurements of each 
configuration and used a jackknife method on these averages to calculate the 
error bars. 

It is clear that the binding energy goes to a non-zero value, with the exception 
of the lowest mass, where the effective binding energy does not seem to converge 
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Figure 6: Like figure ||, but with K2 = 1.259. 



within the hmited distance range. We chose not to display distances larger than 
the position of the maximum of N'{r) in figure 0. 

Unfortunately, the correlation between the mass and the binding energy does 
not appear to be strictly positive. The lowest binding energy belongs to the 
highest mass. We defer a more elaborate discussion to the next section. 

Figure ^ shows the corresponding data for K2 = 1.259 which is the higher of the 
two H2 values used in figure H for N'{r). In this case we used 200 configurations, 
again with 120 origins per configuration. The data look very similar. Here the 
ordering of binding energies follows more clearly that of the constituent masses, 
except for the largest mass mo = 1. 

To see the effect of using smeared sources we show in figure ^ the effective 
binding energy using the unsmeared estimator ( |T7D for k,2 = 1-259, the analogue of 
figure |[ The use of smeared sources (figure || compared to figure 0) does indeed 
improve the convergence to a definite value and the effective binding energies 
are generally smaller. In particular the curve for a bare mass of 1, which keeps 
going down in figure |^ becomes more horizontal. Only the smallest mass is an 
exception, where neither case has converged yet in the distance range shown. 

To determine the binding energies we have fitted the (smeared) effective Eh{r) 
to the form 



Eb{r) = Eb + {c + '-ylnr) r 



(42) 
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Figure 7: Like figure |, but with unsmeared sources. 



which corresponds to the asymptotic forms (JT^) and (|TB|) with c = ln(Z^/Z'^^^) 
and 7 = 2a — /?. In general one would expect c to be positive since 0^ will be 
relatively less effective in creating the two-particle bound state than the single 
particle state. If the single and two-particle propagators behave sufficiently like 
a scalar propagator in flat space we would have a ~ /5 < , hence also 7 < 0. 
The results fo the fits in the region 3 ^ r ^ 20 (recall r^, ~ 21 — 22) are shown 
in the following table. 



mo 


Efe(1.255) 


c(1.255) 


7(1.255) 


^^(1.259) 


c(1.259) 


7(1.259) 


1 

0.316 
0.1 

0.0316 


0.054(1) 
0.078(2) 
0.064(2) 
0.035(2) 


0.053(5) 
0.134(9) 
0.17(1) 
0.17(2) 


-0.100(7) 
-0.16(1) 
-0.09(1) 
+0.03(2) 


0.045(1) 
0.053(1) 
0.038(2) 
0.019(2) 


0.059(4) 
0.19(7) 
0.15(1) 
0.15(2) 


-0.084(6) 
-0.089(9) 
-0.04(1) 
+0.02(2) 



(43) 

The signs of c and 7 appear to follow the above expectations, except for 7 for 
the smallest bare mass. The corresponding effective binding energy is still falling 
for the distances shown and has to be viewed with caution. We cannot improve 
on this because we do not trust distances much beyond the maximum of N'{r) 
since these run into its asymmetric tail. 
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7 Discussion 



The spacetimes produced by the dynamical triangulation model have semiclas- 
sical features near the phase transition. For the binding energy computations 
we have chosen two values of ^2 such that the volume at distance r, N'{r), be- 
haves like that of a four-sphere, for not too large distances r outside a 'planckian 
regime'. 

Let us briefly recall here our scenario for reducing the lattice distance f^: 
increasing we have to tune K2 such that we stay on the same scaling function 
p = rmN'{r)/N. This can be tried for different starting points (/t2,A^), corre- 
sponding to different shapes of p as a function of r/r^. The scaling of p has to 
be carefully re-examined, especially now the phase transition appears to be first 
order. 

The data for a nonzero binding energy are quite convincing and encouraging. 
For given bare mass mo, the resulting renormalized masses m are closer to mo 
and the binding energies turn out to be smaller for the larger k,2 value (1.259), 
than for the smaller value (1.255). This suggests that the renormalized Newton 
constant G is smaller for the larger K2 value, assuming of course G exists and is 
positive for these ^2 values. Since Gq oc 1/^2 and one would indeed expect G to 
decrease when K2 increases. 

We would like to be able to extract G from the data, for example according 
to the nonrelativistic formula 

Eh = 2m-M= ^G^m^ (44) 

This formula is just the familiar energy a^mred/2 of the hydrogen atom in the 
ground state, but with the gravitational parameters substituted as a ^ Grn? 
and the reduced mass mred — »• m/2. Because the formula (|^ ) is nonrelativistic 
it may not suffice to fit the data. To get a rudimentary feeling for corrections to 
( PI) we consider the hamiltonian 

H = 2v/m2 +p2 _ Gm^/r. (45) 

Replacing p — > 1/r and minimising the energy leads to 

Eb = 2m- 2mv/l -G2mV4, (46) 

which suggests that Gm^ = 2 has to be considered 'large'. 

Unfortunately, the data in figure ^ or |^ show no sign of the m^ behaviour of 
(PI), and neither is (|46D of any help. Even the largest constituent mass m ^ 1.2, 
which is evidently so large in lattice units that only qualitative conclusions may 
be drawn from it, leads to a small binding energy. For the lighter constituent 
masses the binding energy is only modestly dependent on m. 
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Perhaps the behaviour for the hghtest constituent mass hints at a possible 
interpretation. If we make the bold assumption that the nonrelativistic formula 
starts making sense for the lightest constituent masses, (|39|,|i3|) lead to a 



renormalised Planck length ^/G = 7.1 and 7.4 respectively for K2 = 1.255 and 
1.259. This appears to contradict the suggestion above that G is smaller for the 
larger K2 value, but the uncertainties in m are greatly magnified by taking its 
fifth power. The values for ^/G are reasonable and they are furthermore similar 
to the (somewhat vague) distance scale the planckian regime ends. On the other 
hand, the size of a nonrelativistic bound state is of order of 1/m, so if indeed 
\/G ~ 7, we should perhaps not be surprised to find odd behaviour for bound 
state sizes of the order of the Planck length, or constituent masses greater than 
the Planck mass 1/VG ~ 0.14. 

Clearly, we are having a problem of separating scales: we would like ^/G <C 
<^ Tm, where = is a measure of the size of our S"^-' universe' 

(the distance where the volume N'{r) is maximal). In our situation, at best, 
\/G ~ m^^ and ^ 21 ^ 3m^^. The size of our universe is only three times 
the Planck length. It is essential for a physical interpretation that as the lattices 
get bigger, the planckian regime and a/G shrink in units of r^- We found some 
evidence for this in our scaling analysis |^ where we ventured a scenario in which 
VG/r^n because of triviality. This approach to zero might then be only 
logarithmic, which would make the problem of scale separation severe from the 
computational point of view. A careful study is needed to clarify these issues. If 
the separation of scales would not materialize we could still study the continuum 
behaviour of the complete function G^'^\r) as a welcome addition to the N'(r) 
used so far. 

Acknowledgements 

We thank Piotr Bialas for useful discussions. The numerical simulations were 
carried out on the IBM SPl at SARA and the Parsytec PowerXplorer at IC3A. 
This work is supported in part by FOM. 

References 

[1] D. Weingarten, Nucl. Phys. B210 [FS6] (1982) 229. 

[2] J. Ambj0rn and J. Jurkiewicz, Phys. Lett. B278 (1992) 42. 

[3] M.E. Agishtein and A.A. Migdal, Mod. Phys. Lett. A7 (1992) 1039. 

[4] G.W. Gibbons, S.W. Hawking and M.J. Perry, Nucl. Phys. B 138 (1978) 
141. 



18 



[5] B.V. de Bakker and J. Smit, Nucl. Phys. B 439 (1995) 239. 

[6] J. Ambj0m and J. Jurkiewicz, Nucl. Phys. B 451 (1995) 643. 

[7] P. Bialas, Z. Burda, A. Krzywicki and B. Petersson, Nucl. Phys. B472 (1996) 
293. 

[8] B.V. de Bakker, Further evidence that the transition of 4D dynamical tri- 
angulation is 1st order, preprint AIAP-1996-040, |hep-lat/9603024i 

[9] J. Frolich, T. Spencer, Commun. Math. Phys. 83 (1982) 411; V. Alessandrini, 
Nucl. Phys. B215 (1983) 337; V. Alessandrini, Ph. Boucaud, Nucl. Phys. 
B225 (1983) 303; K. Jansen, J. Jersak, CB. Lang, T. Neuhaus, G. Vones, 
Nucl. Phys. B265 (1986) 129. 

[10] H.W. Hamber and R.M. Williams, Nucl. Phys. B 435 (1995) 361. 

[11] G. Modanese, Nucl. Phys. B 434 (1995) 697. 

[12] H.W. Hamber and S. Liu, Phys. Lett. B 357 (1995) 51. 

[13] B.V. de Bakker and J. Smit, Nucl. phys. (Proc. Suppl.) 34 (1994) 739. 

[14] B.V. de Bakker and J. Smit, Nucl. Phys. (Proc. Supppl.) 47 (1996) 613. 

[15] S.R. Sharpe, in: CP Violation and the limits of the Standard Model, ed. 
J.F. Donoghue (World Scientific, 1995) p. 377. 

[16] J. Ambj0rn, Z. Burda, J. Jurkiewicz and C.F. Kristjansen, Phys. Rev. D 48 
(1993) 3695. 

[17] B.V. de Bakker and J. Smit, Nucl. Phys. B 454 (1995) 343. 
[18] P. Bialas, private communication. 

[19] I. Antoniadis and E. Mottola, J. Math. Phys. 32 (1991) 1037. 

[20] M.E. Agishtein and A.A. Migdal, Nucl. Phys. B385 (1992) 395. 

[21] J. Ambj0rn, J. Jurkiewicz and C.F. Kristjansen, Nucl. Phys. B393 (1993) 
601. 

[22] B. Briigmann, Phys. Rev. D 47 (1993) 3330. 

[23] S. Catterall, J. Kogut and R. Renken, Phys. Lett. B328 (1994) 277. 
[24] B.V. de Bakker and J. Smit, Phys. Lett. B 334 (1994) 304. 
[25] R. Gupta, Nucl. Phys. B (Proc. Suppl.) 17 (1990) 70. 



19 



